Limitations in harnessing oral RNA interference as an antiviral strategy in Aedes aegypti.

Mosquitoes, particularly Aedes aegypti, are critical vectors for globally significant pathogenic viruses. This study examines the limitations of oral RNA interference (RNAi) as a strategy to disrupt viral transmission by Ae. aegypti. We hypothesized that double-stranded RNA (dsRNA) targeting the Zika virus (ZIKV) or chikungunya virus (CHIKV) genomes produced by engineered bacterial symbionts could trigger an antiviral response. Mosquitoes mono-colonized with Escherichia coli producing dsZIK or dsCHIK did not display reduced viral titers following exposure to virus-contaminated bloodmeals and failed to generate dsZIK- or dsCHIK-derived small interfering RNAs. To address potential limitations of bacterial dsRNA release, we explored dsRNA inoculation via feeding and injection. Although viral replication was impeded in mosquitoes injected with dsZIK or dsCHIK, no antiviral effect was observed in dsRNA-fed mosquitoes. These findings highlight complexities of implementing oral RNAi as an antiviral strategy in Ae. aegypti and warrant further exploration of local and systemic RNAi mechanisms.

Multifaceted contributions of Dicer2 to arbovirus transmission by Aedes aegypti.

Arthropod-borne viruses (arboviruses) transmitted by Aedes aegypti mosquitoes are an increasing threat to global health. The small interfering RNA (siRNA) pathway is considered the main antiviral immune pathway of insects, but its effective impact on arbovirus transmission is surprisingly poorly understood. Here, we use CRISPR-Cas9-mediated gene editing in vivo to mutate Dicer2, a gene encoding the RNA sensor and key component of the siRNA pathway. The loss of Dicer2 enhances early viral replication and systemic viral dissemination of four medically significant arboviruses (chikungunya, Mayaro, dengue, and Zika viruses) representing two viral families. However, Dicer2 mutants and wild-type mosquitoes display overall similar levels of vector competence. In addition, Dicer2 mutants undergo significant virus-induced mortality during infection with chikungunya virus. Together, our results define a multifaceted role for Dicer2 in the transmission of arboviruses by Ae. aegypti mosquitoes and pave the way for further mechanistic investigations.

The E2 glycoprotein holds key residues for Mayaro virus adaptation to the urban Aedes aegypti mosquito.

Adaptation to mosquito vectors suited for transmission in urban settings is a major driver in the emergence of arboviruses. To better anticipate future emergence events, it is crucial to assess their potential to adapt to new vector hosts. In this work, we used two different experimental evolution approaches to study the adaptation process of an emerging alphavirus, Mayaro virus (MAYV), to Ae. aegypti, an urban mosquito vector of many other arboviruses. We identified E2-T179N as a key mutation increasing MAYV replication in insect cells and enhancing transmission after escaping the midgut of live Ae. aegypti. In contrast, this mutation decreased viral replication and binding in human fibroblasts, a primary cellular target of MAYV in humans. We also showed that MAYV E2-T179N generates reduced viremia and displays less severe tissue pathology in vivo in a mouse model. We found evidence in mouse fibroblasts that MAYV E2-T179N is less dependent on the Mxra8 receptor for replication than WT MAYV. Similarly, exogenous expression of human apolipoprotein receptor 2 and Mxra8 enhanced WT MAYV replication compared to MAYV E2-T179N. When this mutation was introduced in the closely related chikungunya virus, which has caused major outbreaks globally in the past two decades, we observed increased replication in both human and insect cells, suggesting E2 position 179 is an important determinant of alphavirus host-adaptation, although in a virus-specific manner. Collectively, these results indicate that adaptation at the T179 residue in MAYV E2 may result in increased vector competence-but coming at the cost of optimal replication in humans-and may represent a first step towards a future emergence event.

Ribosomal RNA (rRNA) sequences from 33 globally distributed mosquito species for improved metagenomics and species identification.

Total RNA sequencing (RNA-seq) is an important tool in the study of mosquitoes and the RNA viruses they vector as it allows assessment of both host and viral RNA in specimens. However, there are two main constraints. First, as with many other species, abundant mosquito ribosomal RNA (rRNA) serves as the predominant template from which sequences are generated, meaning that the desired host and viral templates are sequenced far less. Second, mosquito specimens captured in the field must be correctly identified, in some cases to the sub-species level. Here, we generate mosquito rRNA datasets which will substantially mitigate both of these problems. We describe a strategy to assemble novel rRNA sequences from mosquito specimens and produce an unprecedented dataset of 234 full-length 28S and 18S rRNA sequences of 33 medically important species from countries with known histories of mosquito-borne virus circulation (Cambodia, the Central African Republic, Madagascar, and French Guiana). These sequences will allow both physical and computational removal of rRNA from specimens during RNA-seq protocols. We also assess the utility of rRNA sequences for molecular taxonomy and compare phylogenies constructed using rRNA sequences versus those created using the gold standard for molecular species identification of specimens-the mitochondrial cytochrome c oxidase I (COI) gene. We find that rRNA- and COI-derived phylogenetic trees are incongruent and that 28S and concatenated 28S+18S rRNA phylogenies reflect evolutionary relationships that are more aligned with contemporary mosquito systematics. This significant expansion to the current rRNA reference library for mosquitoes will improve mosquito RNA-seq metagenomics by permitting the optimization of species-specific rRNA depletion protocols for a broader range of species and streamlining species identification by rRNA sequence and phylogenetics.

Innate immune pathways act synergistically to constrain RNA virus evolution in Drosophila melanogaster.

Host-pathogen interactions impose recurrent selective pressures that lead to constant adaptation and counter-adaptation in both competing species. Here, we sought to study this evolutionary arms-race and assessed the impact of the innate immune system on viral population diversity and evolution, using Drosophila melanogaster as model host and its natural pathogen Drosophila C virus (DCV). We isogenized eight fly genotypes generating animals defective for RNAi, Imd and Toll innate immune pathways as well as pathogen-sensing and gut renewal pathways. Wild-type or mutant flies were then orally infected with DCV and the virus was serially passaged ten times via reinfection in naive flies. Viral population diversity was studied after each viral passage by high-throughput sequencing and infection phenotypes were assessed at the beginning and at the end of the evolution experiment. We found that the absence of any of the various immune pathways studied increased viral genetic diversity while attenuating virulence. Strikingly, these effects were observed in a range of host factors described as having mainly antiviral or antibacterial functions. Together, our results indicate that the innate immune system as a whole and not specific antiviral defence pathways in isolation, generally constrains viral diversity and evolution.

APOBEC3C S188I Polymorphism Enhances Context-Specific Editing of Hepatitis B Virus Genome.

Single-nucleotide polymorphism in APOBEC3C (resulting in a serine to isoleucine in position 188) is present in approximately 10% of African populations and greatly enhances restriction against human immunodeficiency virus-1 and simian immunodeficiency virus by improving dimerization and DNA processivity of the enzyme. In this study, we demonstrated in culture and in infected patients that hepatitis B virus (HBV) could be edited by APOBEC3CS188I. Using next-generation sequencing, we demonstrated that APOBEC3CS188I led to enhanced editing activity in 5'TpCpA→5'TpTpA context. This constitutes a new hallmark of this enzyme, which could be used to determine its impact on HBV or nuclear DNA.

Innovative Toolbox for the Quantification of Drosophila C Virus, Drosophila A Virus, and Nora Virus.

Quantification of viral replication underlies investigations into host-virus interactions. In Drosophila melanogaster, persistent infections with Drosophila C virus, Drosophila A virus, and Nora virus are commonly observed in nature and in laboratory fly stocks. However, traditional endpoint dilution assays to quantify infectious titers are not compatible with persistently infecting isolates of these viruses that do not cause cytopathic effects in cell culture. Here we present a novel assay based on immunological detection of Drosophila C virus infection that allows quantification of infectious titers for a wider range of Drosophila C virus isolates. We also describe strand specific RT-qPCR assays for quantification of viral negative strand RNA produced during Drosophila C virus, Drosophila A virus, and Nora virus infection. Finally, we demonstrate the utility of these assays for quantification of viral replication during oral infections and persistent infections with each virus.

Herpes Simplex Virus Type 1 Infection Disturbs the Mitochondrial Network, Leading to Type I Interferon Production through the RNA Polymerase III/RIG-I Pathway.

Viruses have evolved a plethora of mechanisms to impair host innate immune responses. Herpes simplex virus type 1 (HSV-1), a double-stranded linear DNA virus, impairs the mitochondrial network and dynamics predominantly through the UL12.5 gene. We demonstrated that HSV-1 infection induced a remodeling of mitochondrial shape, resulting in a fragmentation of the mitochondria associated with a decrease in their volume and an increase in their sphericity. This damage leads to the release of mitochondrial DNA (mtDNA) to the cytosol. By generating a stable THP-1 cell line expressing the DNase I-mCherry fusion protein and a THP-1 cell line specifically depleted of mtDNA upon ethidium bromide treatment, we showed that cytosolic mtDNA contributes to type I interferon and APOBEC3A upregulation. This was confirmed by using an HSV-1 strain (KOS37 UL98-SPA) with a deletion of the UL12.5 gene that impaired its ability to induce mtDNA stress. Furthermore, by using an inhibitor of RNA polymerase III, we demonstrated that upon HSV-1 infection, cytosolic mtDNA enhanced type I interferon induction through the RNA polymerase III/RIG-I pathway. APOBEC3A was in turn induced by interferon. Deep sequencing analyses of cytosolic mtDNA mutations revealed an APOBEC3A signature predominantly in the 5'TpCpG context. These data demonstrate that upon HSV-1 infection, the mitochondrial network is disrupted, leading to the release of mtDNA and ultimately to its catabolism through APOBEC3-induced mutations. IMPORTANCE Herpes simplex virus 1 (HSV-1) impairs the mitochondrial network through the viral protein UL12.5. This leads to the fusion of mitochondria and simultaneous release of mitochondrial DNA (mtDNA) in a mouse model. We have shown that released mtDNA is recognized as a danger signal, capable of stimulating signaling pathways and inducing the production of proinflammatory cytokines. The expression of the human cytidine deaminase APOBEC3A is highly upregulated by interferon responses. This enzyme catalyzes the deamination of cytidine to uridine in single-stranded DNA substrates, resulting in the catabolism of edited DNA. Using human cell lines deprived of mtDNA and viral strains deficient in UL12, we demonstrated the implication of mtDNA in the production of interferon and APOBEC3A expression during viral infection. We have shown that HSV-1 induces mitochondrial network fragmentation in a human model and confirmed the implication of RNA polymerase III/RIG-I signaling in the capture of cytosolic mtDNA.

Sustained high expression of multiple APOBEC3 cytidine deaminases in systemic lupus erythematosus.

APOBEC3 (A3) enzymes are best known for their role as antiviral restriction factors and as mutagens in cancer. Although four of them, A3A, A3B, A3F and A3G, are induced by type-1-interferon (IFN-I), their role in inflammatory conditions is unknown. We thus investigated the expression of A3, and particularly A3A and A3B because of their ability to edit cellular DNA, in Systemic Lupus Erythematosus (SLE), a chronic inflammatory disease characterized by high IFN-α serum levels. In a cohort of 57 SLE patients, A3A and A3B, but also A3C and A3G, were upregulated ~ 10 to 15-fold (> 1000-fold for A3B) compared to healthy controls, particularly in patients with flares and elevated serum IFN-α levels. Hydroxychloroquine, corticosteroids and immunosuppressive treatment did not reverse A3 levels. The A3AΔ3B polymorphism, which potentiates A3A, was detected in 14.9% of patients and in 10% of controls, and was associated with higher A3A mRNA expression. A3A and A3B mRNA levels, but not A3C or A3G, were correlated positively with dsDNA breaks and negatively with lymphopenia. Exposure of SLE PBMCs to IFN-α in culture induced massive and sustained A3A levels by 4 h and led to massive cell death. Furthermore, the rs2853669 A > G polymorphism in the telomerase reverse transcriptase (TERT) promoter, which disrupts an Ets-TCF-binding site and influences certain cancers, was highly prevalent in SLE patients, possibly contributing to lymphopenia. Taken together, these findings suggest that high baseline A3A and A3B levels may contribute to cell frailty, lymphopenia and to the generation of neoantigens in SLE patients. Targeting A3 expression could be a strategy to reverse cell death and the generation of neoantigens.

Evidence For Long-Lasting Transgenerational Antiviral Immunity in Insects.

Transgenerational immune priming (TGIP) allows memory-like immune responses to be transmitted from parents to offspring in many invertebrates. Despite increasing evidence for TGIP in insects, the mechanisms involved in the transfer of information remain largely unknown. Here, we show that Drosophila melanogaster and Aedes aegypti transmit antiviral immunological memory to their progeny that lasts throughout generations. We observe that TGIP, which is virus and sequence specific but RNAi independent, is initiated by a single exposure to disparate RNA viruses and also by inoculation of a fragment of viral double-stranded RNA. The progeny, which inherit a viral DNA that is only a fragment of the viral RNA used to infect the parents, display enriched expression of genes related to chromatin and DNA binding. These findings represent a demonstration of TGIP for RNA viruses in invertebrates, broadly increasing our understanding of the immune response, host genome plasticity, and antiviral memory of the germline.

Non-retroviral Endogenous Viral Element Limits Cognate Virus Replication in Aedes aegypti Ovaries.

Endogenous viral elements (EVEs) are viral sequences integrated in host genomes. A large number of non-retroviral EVEs was recently detected in Aedes mosquito genomes, leading to the hypothesis that mosquito EVEs may control exogenous infections by closely related viruses. Here, we experimentally investigated the role of an EVE naturally found in Aedes aegypti populations and derived from the widespread insect-specific virus, cell-fusing agent virus (CFAV). Using CRISPR-Cas9 genome editing, we created an Ae. aegypti line lacking the CFAV EVE. Absence of the EVE resulted in increased CFAV replication in ovaries, possibly modulating vertical transmission of the virus. Viral replication was controlled by targeting of viral RNA by EVE-derived P-element-induced wimpy testis-interacting RNAs (piRNAs). Our results provide evidence that antiviral piRNAs are produced in the presence of a naturally occurring EVE and its cognate virus, demonstrating a functional link between non-retroviral EVEs and antiviral immunity in a natural insect-virus interaction.

Differential Small RNA Responses against Co-Infecting Insect-Specific Viruses in Aedes albopictus Mosquitoes.

The mosquito antiviral response has mainly been studied in the context of arthropod-borne virus (arbovirus) infection in female mosquitoes. However, in nature, both female and male mosquitoes are frequently infected with insect-specific viruses (ISVs). ISVs are capable of infecting the reproductive organs of both sexes and are primarily maintained by vertical transmission. Since the RNA interference (RNAi)-mediated antiviral response plays an important antiviral role in mosquitoes, ISVs constitute a relevant model to study sex-dependent antiviral responses. Using a naturally generated viral stock containing three distinct ISVs, Aedes flavivirus (AEFV), Menghai rhabdovirus (MERV), and Shinobi tetra virus (SHTV), we infected adult Aedes albopictus females and males and generated small RNA libraries from ovaries, testes, and the remainder of the body. Overall, both female and male mosquitoes showed unique small RNA profiles to each co-infecting ISV regardless of the sex or tissue tested. While all three ISVs generated virus-derived siRNAs, only MERV generated virus-derived piRNAs. We also studied the expression of PIWI genes in reproductive tissues and carcasses. In contrast to Piwi5-9, Piwi1-4 were abundantly expressed in ovaries and testes, suggesting that Piwi5-9 are involved in exogenous viral piRNA production. Together, our results show that ISV-infected Aedes albopictus produce viral small RNAs in a virus-specific manner and that male mosquitoes mount a similar small RNA-mediated antiviral response to that of females.

Tudor-SN Promotes Early Replication of Dengue Virus in the Aedes aegypti Midgut.

Diseases caused by mosquito-borne viruses have been on the rise for the last decades, and novel methods aiming to use laboratory-engineered mosquitoes that are incapable of carrying viruses have been developed to reduce pathogen transmission. This has stimulated efforts to identify optimal target genes that are naturally involved in mosquito antiviral defenses or required for viral replication. Here, we investigated the role of a member of the Tudor protein family, Tudor-SN, upon dengue virus infection in the mosquito Aedes aegypti. Tudor-SN knockdown reduced dengue virus replication in the midgut of Ae. aegypti females. In immunofluorescence assays, Tudor-SN localized to the nucleolus in both Ae. aegypti and Aedes albopictus cells. A reporter assay and small RNA profiling demonstrated that Tudor-SN was not required for RNA interference function in vivo. Collectively, these results defined a novel proviral role for Tudor-SN upon early dengue virus infection of the Ae. aegypti midgut.

Host nutritional status affects alphavirus virulence, transmission, and evolution.

Malnourishment, specifically overweight/obesity and undernourishment, affects more than 2.5 billion people worldwide, with the number affected ever-increasing. Concurrently, emerging viral diseases, particularly those that are mosquito-borne, have spread dramatically in the past several decades, culminating in outbreaks of several viruses worldwide. Both forms of malnourishment are known to lead to an aberrant immune response, which can worsen disease outcomes and reduce vaccination efficacy for viral pathogens such as influenza and measles. Given the increasing rates of malnutrition and spread of arthropod-borne viruses (arboviruses), there is an urgent need to understand the role of host nutrition on the infection, virulence, and transmission of these viruses. To address this gap in knowledge, we infected lean, obese, and undernourished mice with arthritogenic arboviruses from the genus Alphavirus and assessed morbidity, virus replication, transmission, and evolution. Obesity and undernourishment did not consistently influence virus replication in the blood of infected animals except for reductions in virus in obese mice late in infection. However, morbidity was increased in obese mice under all conditions. Using Mayaro virus (MAYV) as a model arthritogenic alphavirus, we determined that both obese and undernourished mice transmit virus less efficiently to mosquitoes than control (lean) mice. In addition, viral genetic diversity and replicative fitness were reduced in virus isolated from obese compared to lean controls. Taken together, nutrition appears to alter the course of alphavirus infection and should be considered as a critical environmental factor during outbreaks.

Chikungunya Virus Vaccine Candidates with Decreased Mutational Robustness Are Attenuated In Vivo and Have Compromised Transmissibility.

Chikungunya virus (CHIKV) is a reemerged arbovirus, a member of the Togaviridae family. It circulates through mosquito vectors mainly of the Aedes family and a mammalian host. CHIKV causes chikungunya fever, a mild to severe disease characterized by arthralgia, with some fatal outcomes described. In the past years, several outbreaks mainly caused by enhanced adaptation of the virus to the vector and ineffective control of the contacts between infected mosquito populations and the human host have been reported. Vaccines represent the best solution for the control of insect-borne viruses, including CHIKV, but are often unavailable. We designed live attenuated CHIKVs by applying a rational genomic design based on multiple replacements of synonymous codons. In doing so, the virus mutational robustness (capacity to maintain phenotype despite introduction of mutations to genotype) is decreased, driving the viral population toward deleterious evolutionary trajectories. When the candidate viruses were tested in the insect and mammalian hosts, we observed overall strong attenuation in both and greatly diminished signs of disease. Moreover, we found that the vaccine candidates elicited protective immunity related to the production of neutralizing antibodies after a single dose. During an experimental transmission cycle between mosquitoes and naive mice, vaccine candidates could be transmitted by mosquito bite, leading to asymptomatic infection in mice with compromised dissemination. Using deep-sequencing technology, we observed an increase in detrimental (stop) codons, which confirmed the effectiveness of this genomic design. Because the approach involves hundreds of synonymous modifications to the genome, the reversion risk is significantly reduced, rendering the viruses promising vaccine candidates.IMPORTANCE Chikungunya fever is a debilitating disease that causes severe pain to the joints, which can compromise the patient's lifestyle for several months and even in some grave cases lead to death. The etiological agent is chikungunya virus, an alphavirus transmitted by mosquito bite. Currently, there are no approved vaccines or treatments against the disease. In our research, we developed novel live attenuated vaccine candidates against chikungunya virus by applying an innovative genomic design. When tested in the insect and mammalian host, the vaccine candidates did not cause disease, elicited strong protection against further infection, and had low risk of reversion to pathogenic phenotypes.

Seasonal Genetic Drift of Human Influenza A Virus Quasispecies Revealed by Deep Sequencing.

After a pandemic wave in 2009 following their introduction in the human population, the H1N1pdm09 viruses replaced the previously circulating, pre-pandemic H1N1 virus and, along with H3N2 viruses, are now responsible for the seasonal influenza type A epidemics. So far, the evolutionary potential of influenza viruses has been mainly documented by consensus sequencing data. However, like other RNA viruses, influenza A viruses exist as a population of diverse, albeit related, viruses, or quasispecies. Interest in this quasispecies nature has increased with the development of next generation sequencing (NGS) technologies that allow a more in-depth study of the genetic variability. NGS deep sequencing methodologies were applied to determine the whole genome genetic heterogeneity of the three categories of influenza A viruses that circulated in humans between 2007 and 2012 in France, directly from clinical respiratory specimens. Mutation frequencies and single nucleotide polymorphisms were used for comparisons to address the level of natural intrinsic heterogeneity of influenza A viruses. Clear differences in single nucleotide polymorphism profiles between seasons for a given subtype also revealed the constant genetic drift that human influenza A virus quasispecies undergo.

Immune priming and clearance of orally acquired RNA viruses in Drosophila.

Immune responses in insects are differentially triggered depending on the infection route used by the pathogen. In most studies involving Drosophila melanogaster and viruses, infection is done by injection, while oral infection, which is probably the most common route of viral entry in nature, remains unexplored. Here, we orally infected adults and larvae from wild-type and RNA interference (RNAi) mutant flies with different RNA viruses. We found that, in contrast with what is observed following virus injection, oral infections initiated at larval or adult stages are cleared in adult flies. Virus elimination occurred despite a larger infectious dose than for injected flies and evidence of viral replication. RNAi mutant flies suffered greater mortality relative to wild-type flies following oral infection, but they also eliminated the virus, implying that RNAi is not essential for viral clearance and that other immune mechanisms act during oral infections. We further showed that information of infection by RNA viruses acquired orally leaves a trace under a DNA form, which confers protection against future reinfection by the same virus. Together, this work presents evidence of clearance and immune priming for RNA viruses in insects and challenges the current view of antiviral immunity in insects.

Dicer-2-Dependent Generation of Viral DNA from Defective Genomes of RNA Viruses Modulates Antiviral Immunity in Insects.

The RNAi pathway confers antiviral immunity in insects. Virus-specific siRNA responses are amplified via the reverse transcription of viral RNA to viral DNA (vDNA). The nature, biogenesis, and regulation of vDNA are unclear. We find that vDNA produced during RNA virus infection of Drosophila and mosquitoes is present in both linear and circular forms. Circular vDNA (cvDNA) is sufficient to produce siRNAs that confer partially protective immunity when challenged with a cognate virus. cvDNAs bear homology to defective viral genomes (DVGs), and DVGs serve as templates for vDNA and cvDNA synthesis. Accordingly, DVGs promote the amplification of vDNA-mediated antiviral RNAi responses in infected Drosophila. Furthermore, vDNA synthesis is regulated by the DExD/H helicase domain of Dicer-2 in a mechanism distinct from its role in siRNA generation. We suggest that, analogous to mammalian RIG-I-like receptors, Dicer-2 functions like a pattern recognition receptor for DVGs to modulate antiviral immunity in insects.

Attenuation of RNA viruses by redirecting their evolution in sequence space.

RNA viruses pose serious threats to human health. Their success relies on their capacity to generate genetic variability and, consequently, on their adaptive potential. We describe a strategy to attenuate RNA viruses by altering their evolutionary potential. We rationally altered the genomes of Coxsackie B3 and influenza A viruses to redirect their evolutionary trajectories towards detrimental regions in sequence space. Specifically, viral genomes were engineered to harbour more serine and leucine codons with nonsense mutation targets: codons that could generate Stop mutations after a single nucleotide substitution. Indeed, these viruses generated more Stop mutations both in vitro and in vivo, accompanied by significant losses in viral fitness. In vivo, the viruses were attenuated, generated high levels of neutralizing antibodies and protected against lethal challenge. Our study demonstrates that cornering viruses in 'risky' areas of sequence space may be implemented as a broad-spectrum vaccine strategy against RNA viruses.

Uncovering the Repertoire of Endogenous Flaviviral Elements in Aedes Mosquito Genomes.

Endogenous viral elements derived from nonretroviral RNA viruses have been described in various animal genomes. Whether they have a biological function, such as host immune protection against related viruses, is a field of intense study. Here, we investigated the repertoire of endogenous flaviviral elements (EFVEs) in Aedes mosquitoes, the vectors of arboviruses such as dengue and chikungunya viruses. Previous studies identified three EFVEs from Aedes albopictus cell lines and one from Aedes aegypti cell lines. However, an in-depth characterization of EFVEs in wild-type mosquito populations and individual mosquitoes in vivo has not been performed. We detected the full-length DNA sequence of the previously described EFVEs and their respective transcripts in several A. albopictus and A. aegypti populations from geographically distinct areas. However, EFVE-derived proteins were not detected by mass spectrometry. Using deep sequencing, we detected the production of PIWI-interacting RNA-like small RNAs, in an antisense orientation, targeting the EFVEs and their flanking regions in vivo The EFVEs were integrated in repetitive regions of the mosquito genomes, and their flanking sequences varied among mosquito populations. We bioinformatically predicted several new EFVEs from a Vietnamese A. albopictus population and observed variation in the occurrence of those elements among mosquitoes. Phylogenetic analysis of an A. aegypti EFVE suggested that it integrated prior to the global expansion of the species and subsequently diverged among and within populations. The findings of this study together reveal the substantial structural and nucleotide diversity of flaviviral integrations in Aedes genomes. Unraveling this diversity will help to elucidate the potential biological function of these EFVEs.IMPORTANCE Endogenous viral elements (EVEs) are whole or partial viral sequences integrated in host genomes. Interestingly, some EVEs have important functions for host fitness and antiviral defense. Because mosquitoes also have EVEs in their genomes, characterizing these EVEs is a prerequisite for their potential use to manipulate the mosquito antiviral response. In the study described here, we focused on EVEs related to the Flavivirus genus, to which dengue and Zika viruses belong, in individual Aedes mosquitoes from geographically distinct areas. We show the existence in vivo of flaviviral EVEs previously identified in mosquito cell lines, and we detected new ones. We show that EVEs have evolved differently in each mosquito population. They produce transcripts and small RNAs but not proteins, suggesting a function at the RNA level. Our study uncovers the diverse repertoire of flaviviral EVEs in Aedes mosquito populations and contributes to an understanding of their role in the host antiviral system.